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A review  of  all  the  currently  available  papers  concerned  with  the  solution 

of  the  problem  of  the  induced  velocity  field  of  a helicopter  rotor  is  presented. 

o 

The  assumptions,  applications,  and  limitations  of  each  paper  are  given.  Since, 
in  general,  all  the  theories  are  very  limited  in  their  applicability,  it  is 
recommended  that  further  theoretical  and  experimental  investigations  be 
initiated  which  will  be  directed  toward  a more  accurate  prediction  of  the 
induced  velocity  field  about  a rotor. 
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Introduction 

Although  the  concept  of  the  lifting  rotor  is  quite  old,  it  has  not  been 
until  recently  that  the  special  problems  involving  the  air  flow  through  and 
around  such  a rotor  have  been  theoretically  and  experimentally  investigated.  In 
fact,  in  general,  such  investigations  were  necessarily  delayed  until  the 
successful  flight  of  a lifting  rotor  in  various  regimes  began  to  point  out 
control  and  vibrational  difficulties  which  were  obviously  brought  about  by  the 
variation  in  the  rotor  inflow  velocities.  The  theories  available  were  merely 
extensions  of  the  propeller  theories  whose  origin,  except  for  the  vortex 
concept,  lay  in  the  older  marine  screw  theories.  While  these  momentum,  blade 
element,  and  vortex  theories  are  adequate  for  performance  estimation,  they  are 
woefully  inadequate  for  use  in  predicting  blade  motions  for  stability  and 
control,  vibration,  and  blade  stress  problems.  This  is  particularly  true  at 
slow  forward  flight  speeds  and  vertical  descent. 

The  purpose  of  this  review  then  is  to  trace  the  development  of  the  rotor 
induced  velocity  field  theory,  to  present  the  applicability  and  limitations  of 
each  theory,  to  reach  a conclusion  concerning  the  state  of  the  knowledge  and  to 
point  out,  if  possible,  ways  and  means  of  improving  or  extending  the  theory. 

This  review  is  divided  into  four  main  parts  as  fellows  % momentum  theory, 
blade  element  theory,  combined  momentum  and  blade  element  theory,  and  vortex 
theory.  The  development  of  each  concept  is  treated  chronologically  from  its 
beginning  to  its  present  state.  Since  the  investigations  and  theoretical 
analyses  of  marine  screws,  propellers,  and  rotors  prior  to  1935  are  adequately 
presented  and  referenced  by  H.  Glauert  in  reference  1,  this  work  is  the  primary 
source  for  material  prior  to  that  year..  Thereafter,  only  the  individual  papers 
that  dealt  solely  with  rotors  and  their  special  problems  are  reviewed. 
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This  literature  review  of  rotor  induced  velocity  field  theories  was 
carried  out  at  Princeton  University  under  the  sponsorship  and  with  the  financial 
assistance  of  the  Office  of  Naval  Research.  ^ 

It  night  be  well  to  briefly  review  the  various  flight  regimes  of  the  '] 

helicopter  before  starting  the  review.  There  are  three  basic  flow  patterns! 

ii 

The  first  corresponds  to  the  propeller  working  state  and  includes  hovering, 
horizontal,  and  ascending  flight.  In  these  flight  regimes,  the  air  at  a great 
distance  ahead  of  the  rotor  and  the  air  passion  through  the  rotpr  in  general 
proceed  in  a downward  direction  with  respect  to  the  tip-path  plane.  The  second 
flow  pattern  is  designated  the  windmill -brake  state.  This  corresponds  to  a 
flight  condition  in  which  the  rotor  is  disconnected  from  the  engine  and  is 
driven  by  the  air  passing  through  it.  In  this  case,  the  air  at  a great  distance 
ahead  of  the  rotor  and  the  air  passing  through  the  rotor  proceed  in  an  upward 
direction  with  respect  to  the  tip  path  plane.  The  third  flow  pattern  lies 
between  the  first  two.  This  corresponds  to  certain  portions  of  the  descending 
flight  regime  and  is  called  the  vortex-ring  state.  Rotor  operation  in  this 
state  is  usually  very  rough.  The  air  at  a great  distance  ahead  of  the  rotor  is 
moving  upward  while  the  air  passing  through  the  rotor  disk  is  moving  downward 
with  respect  to  the  tip  path  plane.  The  limiting  condition  between  the  vortex- 
ring and  the  windmill-brake  states  is  called  the  ideal  autorotation  condition 
in  which  very  little  air  passes  through  the  disk.  The  limit  between  the 
propeller-working  and  vortex-ring  states  in  vertical  flight  is  of  course  the 
hovering  condition.  In  forward  flight  this  limit  will  lie  along  some  descending 
flight  path  since  the  tip  path  plane  is  usually  inclined  forward. 
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List  of  Symbols 


Throughout  this  paper  an  effort  was  made  to  use  as  many  as  possible  of  the 
standard  symbols  for  helicopters  as  set  forth  in  reference  2.  This  necessitated 
the  replacement  of  some  of  the  symbols  in  some  of  the  papers  reviewed  with  these 
standard  symbols.  In  sane  cases,  however,  this  was  not  feasible  and  some 
duplication  exists.  In  these  cases,  the  symbol  is  defined  where  used,  and  it 
is  hoped  that  no  confusion  will  arise. 


slope  of  lift  curve  for  blade  element,  per  radian  (blade  element  section) 
axial  interference  factor,  = !/"  (1  + a)  (momentum  section) 

rotational  interference  factor,  u*  * Z&'fl 

longitudinal  angle  between  the  rotor  shaft  and  the  perpendicular  to  the 
tip  path  plane,  radians 

number  of  blades  in  rotor 

blade  chord  at  radlu6  r 

section  profile  - drag  coefficient 

section  lift  coefficient 


rotor  thrust  coefficient 


ViUUDW  LWlllWACUU  if  ll*  J?** 

rotor  thrust  coefficient  based  on  flight  velocity  — y2 
non-dimensional  function  f (<P,  b ) in  Goldstein  vortex  theory 


lift,  pounds  . 

static  pressure,  pounds  per  square  foot 
power,  foot-pounds  per  second 
rotor  torque,  foot-pounds 

radius  of  blade  element  from  rotor  center,  feet 

radius  of  rotor  blade  tip,  feet 

radius  of  ultimate  rotor  wake,  feet 

gap  distance  between  vortex  sheets  forming  wake,  feet 


rotor  thrust,  pounds 

rotational  component  of  wake  velocity,  radians  per  second 

component  of  resultant  velocity  at  blade  element  that  Is  normal  to  blade 
axis,  feet  per  second 

component  of  resultant  velocity  at  blade  element  that  Is  normal  to  blade 
span  axis  and  Uj.  , feet  per  second 

coiqponent  of  resultant  velocity  in  ultimate  wake  that  is  normal  to  a 
plane  parallel  to  tip  path  plane,  feet  per  second 

component  of  resultant  velocity  at  blade  element  that  is  normal  to  blade 
span  axis  and  to  axis  of  no-feathering,  feet  per  second 

normal  component  of  induced  velocity  at  tip-path  plane,  feet  per  second 

vector  sum  of  the  three  components  of  induced  velocity,  feet  per  second 

average  value  of  fore-and-aft  induced  velocity  at  the  center  of  the  disk 
feet  per  second 

rate  of  change  of  induced  velocity  at  center  of  disk,  feet  per  second 

radial  component  of  velocity  induced  at  a point  P by  a vortex  ring,  feet 
per  second 

axial  component  of  velocity  induced  at  a point  P by  a vortex  ring,  feet 
per  second 

axial  component  of  induced  velocity  at  3/h  radius  due  to  bound 
vortices,  feet  per  second 

axial  component  of  induced  velocity  at  rotor  center  due  to  ro  - 
component  of  the  free  tip  vortices,  feet  per  second 

axial  component  of  induced  velocity  at  3/h  radius  due  to  ro 
component  of  the  free  tip  vortices,  feet  per  second 

axial  component  of  induced  velocity  at  3/U  radius  due  to  the  free- 
radial  and  tip  vortices,  feet  per  second 

true  airspeed  of  helicopter  along  flight  path,  feet  per  second 

component  of  true  airspeed  of  helicopter  normal  to  tip  path  plane,  feet 
per  second 

rotational  component  of  induced  velocity,  feet  per  second 
gross  weight  of  helicopter,  pounds 

velocity  of  axial  translation  of  the  rigid  helicoidal  surface,  feet  per 
second 
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x ratio  of  blade-element  radius  to  rotor-blade  radius,  r/R 

oc,  attitude  angle  of  fuselage;  angle  between  rotor  axis  and  the  vertical, 
radians 

oCy.  blade-element  single  of  attack  measured  from  line  of  zero  lift,  radians 

ccs  angle  between  the  slipstream  and  the  perpendicular  to  the  tip-path  plane 
at  the  rotor  disk,  radians 

ccT  angle  of  attack  of  tip-path  plane  measured  in  the  longitudinal  plane 
between  the  flight-path  velocity  vector  and  tip-path  plane,  radians 

P circulation  of  blade  element  at  radius  r and  azimuth  angle 

constants  in  expression  for  T (r  - j/v  ir) 

Cj 

£ ratio  of  profile-drag  coefficient  to  lift  coefficient  -j— 

■>j  propeller  efficiency 

© blade  pitch  angle  at  particular  blade  radius  and  azimuth  position,  radians 

j <C  coefficient  in  Goldstein  vortex  theory  as  modified  by  Lock  which  represents 

the  circulation  loss  along  the  blade 

j*.y  in-plane  velocity  ratio  at  tip  path  plane 

f>  density  of  air,  slugs  per  cubic  foot 

<TJl  rotor  solidity  at  radius,  r 

<p  inflow  angle  at  blade  element  measured  in  plane  perpendicular  to  blade 
axis  and  between  tip-path  plane  and  relative  wind,  radians 

<p  angle  that  the  helical  surfaces  of  the  vortex  sheets  make  with  the 
' propeller  disk  plane,  radians 

potential  at  a point  P due  to  a closed  vortex  ring 

X wake  skew  angle,  angle  between  perpendicular  to  tip  path  plane  and  wake 
boundary  at  a far  distance  below  the  rotor  disk  in  a longitudinal  plane, 
radians 

~'fr  azimuth  position  of  blade  measured  from  downwind  position  in  direction  of 
rotation,  radians 

cu  solid  angle  subtended  at  a point  by  a vortex  ring 
/l  rotor  angular  velocity,  radians  per  second 

Subscript  - r,  unless  otherwise  noted,  denotes  radial  position  on  blade 
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Momentum  Theory 


The  simple  axial  momentum  theory  as  first  originated  by  W.  J.  M.  Rankine 
in  1865  (reference  3)  was  based  on  the  concept  that  the  forces  experienced  by 
a set  of  rotating  blades  immersed  in  a frictionless,  incompressible  fluid  were 
equal  to  that  necessary  to  provide  the  resulting  motions  of  the  fluid.  In 
further  simplifying  the  physical  picture,  it  was  assumed  that  the  slipstream 
had  no  rotational  component,  that  the  thrust  of  the  propeller  was  constant 
over  the  disk,  that  the  axial  velocity  of  the  fluid  had  a constant  value  over 
the  disk  and  over  a cross-section  of  the  ultimate  wake,  and  that  the  resultant 
pressure  force  on  the  whole  fluid  was  zero.  Then  assuming  the  propeller  to  be 
at  rest  and  the  fluid  to  have  a uniform  velocity  along  its  axis,  the  thrust  on 
the  propeller  was  equated  to  the  axial  increase  in  momentum  in  unit  time,  so 


T - tt  VY) 

where  the  subscript  1 refers  to  the  ultimate  wake.  Since  the  propeller  was 
assumed  to  be  at  rest,  no  useful  work  was  done  and  the  power  absorbed  was  taken 
as  being  equal  to  the  increase  in  kinetic  energy  of  the  slipstream  in  unit  time, 


so  that 


AK.e.--  P » 


Then  using  the  condition  of  continuity  of  flow,  it  was  shown  that 

Vv  = i(Up(  + Vv) 

This  was  one  of  the  most  important  features  of  the  momentum  theory.  It 
pointed  out  that  the  axial  velocity  of  the  fluid  moving  through  the  propeller 
disk  was  greater  than  the  forward  velocity  of  the  disk.  It  further  showed 
that  the  induced  velocity  at  the  disk  was  one-half  its  value  in  the  ultimate 


II 


I. 

t ■ 


8. 


8. 


wake.  This  result  may  also  be  obtained  by  a suitable  application  of  Bernoulli 5 s 
equation,  the  continuity  equation,  and  the.  axial  momentum  equation.  The 
efficiency  of  propulsion  was  then  shown  to  be 
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By  introducing  an  axial  interference  factor  a,  where 

% sr  \T  ( , + a) 

it  was  shown  that  the  ideal  propeller  efficiency  became 

^ = (I  + a) 

which  is  the  highest  efficiency  that  can  be  obtained  from  a propeller  of  given 
disk  area,  velocity  of  advance,  and  power  absorption . A consequence  of  this 
relation  was  that  the  minimum  loss  of  energy  for  a given  thrust  occurred  when 
the  thrust  was  uniformly  distributed  over  the  whole  disk.  This  distribution 
maintained  a constant  value  of  the  axial  interference  factor. 

R.  E.  Froude,  in  1889  (reference  4-)  introduced  the  actuator  disk  concept. 
This  concept  represented  the  propeller  as  a disk  at  which  there  is  a sudden 
increase  in  pressure  without  a discontinuity  in  velocity.  This  may  be  physically 
represented  by  considering  a close  pair  of  contra-rotating  coaxial  propellers 
having  an  infinite  number  of  frictionless  blades  so  designed  that  the 
rotational  velocity  component  of  the  front  propeller  is  exactly  cancelled  by 
the  component  of  the  rear  for  each  radial  annulus,  and  that  the  blade  angles  are 
chosen  in  such  a manner  as  to  yield  a uniform  distribution  of  thrust  over  the 
disk. 

N.  E.  Joukowski,  in  1918  (reference  5),  and  A.  Betz,  in  1920  (reference  6), 
extended  the  axial  momentum  theory  to  the  more  general  case  by  taking  the 
rotational  velocity  component  into  account.  A similar  derivation  is  given  by 

O 

H.  Glauert,  (reference  l).  It  was  assumed  that  the  actuator  disk  introduced  a 
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rotational  component  to  the  fluid  velocity  while  the  axial  and  radial  component 
remained  unchanged.  This  further  generalization  of  the  theory  resulted  in 
equations  which  were  rather  difficult  to  solve,  unless  the  wake  angular  velocity 
component  was  known  as  a function  of  the  radius.  Thus,  an  exact  solution  could 
be  realized  when  the  flow  in  the  slipstream  was  assumed  to  be  irrotational 
except  along  the  axis.  This  assumption  implied  that  the  rotational  momentum 
had  the  same  value  for  all  radial  elements  and  further  that  the  circulation 
along  the  blade  was  constant.  Such  a distribution  was  not  physically  realizable 
since  it  would  mean  that  close  to  the  axis  the  wake  had  a greater  angular 
velocity  than  the  blade  itself.  Consequently  the  blade  was  usually  assumed  to 
begin  at  the  radial  station  at  which  the  angular  rotation  imported  to  the  slip- 
stream was  identical  to  the  rotational  velocity  of  the  blade  itself.  For  this 
particular  case  it  was  shown  that  when  rotation  of  the  slipstream  was  taken  into 
account,  the  induced  velocity  at  the  disk  was  not  necessarily  one-half  its 
value  in  the  ultimate  wake.  In  fact 

In  general  since  the  angular  velocity  imported  to  the  slipstream  was  very 
much  smaller  than  the  angular  velocity  of  the  propeller  blades,  it  was  possible 
to  simplify  the  general  equations  by  neglecting  certain  terms  involving  the 
square  of  the  wake  rotational  component.  This  yielded  results  similar  to  the 
simple  axial  momentum  theory  for  the  relationship  connecting  the  thrust  and 
axial  velocity.  The  equations  for  the  thrust  and  torque  of  an  ideal  frictionless 
propeller  are  given  on  page  196  of  reference  1 and  are  as  follows: 

AT  = 4tt ( I + aj  a r Jr 
JlQ  = HTt  fVrn(i  +a  )a'r3Jr 

where  a,  the  axial  interference  factor,  is  defined  by 

Vf  = V + a) 
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and.  a',  the  rotational  Interference  factor.  Is  defined  by 

u'  = aa'fl 

where  u‘  is  the  wake  rotational  component.  The  factors  a and  a'  are  related  by 
the  expression 

VY  (i  ta.'ja.  = ( > - A1) 

The  efficiency  for  this  case  is 


The  condition  for  minimum  loss  of  energy  when  the  slipstream  rotation  vas 
included,  neglecting  second  order  and  higher  small  quantities,  vas  for  a 
distribution  of  thrust  such  that  the  efficiency  had  the  same  value  for  all 
elements.  This  simple  condition  applied  only  to  lightly  loaded  propellers. 

The  momentum  theory  may  be  extended  to  include  frictional  drag  of  the 
blades  and  the  interference  of  nearby  bodies.  Rankine,  in  reference  3,  gave  an 
estimate  of  these  effects. 

The  great  disadvantage  of  the  momentum  theory  vas  that  it  gave  no  indication 
of  the  manner  in  which  the  blades  should  be  designed,  other  than  that  the 
diameter  of  the  propeller  should  be  as  large  as  possible. 

A classical  approach  to  the  problem  has  been  given  by  Mangier  in  references 
7,  8,  and  9*  He  also  used  the  actuator  disk  concept  in  which  the  rotor  was 
replaced  by  a circular  disk  with  a pressure  step  across  it.  This  implied  a 
rotor  having  an  infinite  number  of  blades  but  vas  a good  assumption  if  the 
rotor  blades  were  driven  at  a high  angular  velocity.  A second  assumption  vas 
made  in  which  it  vas  assumed  that  the  rotor  vas  either  lightly  loaded  or  that 
the  velocity  along  the  flight  path  vas  large  so  that  the  induced  velocities  were 
small  compared  to  the  flight  path  velocity.  It  vas  also  assumed  that  the 
induced  velocities  that  were  associated  with  the  rotor  torque  could  be  neglected. 
The  second  assumption  allowed  the  author  to  use  the  linearized  theory  as  explained 


»H‘I  Wl  cmif) 


■iBl1litfc.ll,  I Utt-  -j^ 


11. 


by  Burgers  in  reference  10.  The  Euler  Equations  could  then  be  simplified  bo  that 
only  linear  terms  were  retained: 

h 

I 


-vf ?■-*  — p 


where  £ is  the  vector  sum  of  the  three  components  of  induced  velocity  and  p is 
the  static  pressure.  Since  the  continuity  equation  had  to  be  satisfied,  i.e. 

div  v = O 


it  was  shown  that 


div  68ao  Is  - P*/?  = O 


The  static  pressure  function  was  therefore  a potential  function  and  the 

J 

acceleration,  - V-.  was  the  gradient  of  this  potential  function.  Then 

Jx 

by  integrating  the  first  equation  above  it  was  shown  that 

x _ / rx  . 

V “ J S*A0  />  *X 

I'-oO 

The  lower  limit  was  taken  as  +«0  since  the  induced  velocities  must  vanish 


at  a far  distance  ahead  of  the  rotor  and  the  integration  was  performed  for 
constant  y and  * . This  result  applied  everywhere  but  inside  the  wake  since 

the  integration  could  not  be  performed  across  the  disk.  However,  it  was 
explained  that  since  the  integrated  equation  was  everywhere  continuous  outside 
the  disk  and  wake,  the  analytical  continuation  of  this  expression  could  be  used 
inside  the  wake  if  another  function  £ (V,  <)  were  added  to  it  to  make  the 
velocity  continuous  across  the  disk.  Thus  for  the  space  within  the  wake 
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The  wake  was  not,  however,  irrotational  and  a velocity  discontinuity  existed 
at  the  wake  boundaries. 

Mangier* s solution  of  the  Laplace  equation 
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followed  that  of  Kinner  in  reference  11.  This  latter  paper  was  originally 
intended  as  a contribution  to  the  theory  of  the  autdgiro,  but,  in  order  to  limit 
its  scope  and  due  to  mathematical  difficulties,  was  restricted  to  the  problem 
of  the  solid  flat  circular  disk.  Mangier  in  a manner  similar  to  Kinner  solved 
the  Laplace  equation  in  terms  of  Legendre  functions  of  the  elliptic  coordinates 
of  the  rotor  disk  (i.e.  ellipsoidal  harmonics).  These  functions  were  dis- 
continuous between  the  two  faces  of  the  disk  but  were  continuous  everywhere 
else.  It  was  then  possible  to  solve  the  induced  velocity  integral  equation 
but  this  was  difficult  and  lengthy  even  for  simple  distributions.  The  solution 
was  then  restricted  to  axially  symmetric  pressure  distributions  in  which  the 
load  was  a function  of  disk  radius  only  and  not  blade  azimuth  position.  Three 
pressure  distributions  were  used.  The  first  term  of  the  series  of  Legendre 
functions  gave  an  elliptical  distribution  and  the  thrust.  The  second  term 
gave  a moment  but  this  term  was  not  used.  The  third  term  when  combined 
linearly  with  the  first  term  in  a suitable  manner  gave  a pressure  distribution 
which  was  zero  at  the  center  and  circumference  of  the  rotor  disk  and  was  a 
reasonable  approximation  for  the  load  existing  on  a rotor  in  forward  flight. 

The  first  and  combined  first  and  third  pressure  distributions  were  the  ones 
used  in  the  report  for  the  calculation  of  the  induced  velocity  perpendicular 
to  the  disk  and  the  vertical  induced  velocity  at  a point  far  behind  the  disk. 
These  calculations  were  tabulated  and  plotted  for  rotor  angles  of  incidence 
of  0,  15,  30,  45,  and  90  degrees. 

The  results  of  these  calculations  showed  that  the  induced  velocity 
distribution  over  a rotor  disk  was  feu?  from  constant  or  a linear  variation. 
Within  the  limitations  of  the  theory  and  for  the  angles  of  incidence  within 
the  realm  of  possibility  for  a helicopter  in  forward  flight,  the  induced 
velocity  for  the  first  pressure  distributions  was  directed  upward  at  the  front 


and  sides  and  downward  over  the  remainder  of  the  disk.  For  the  combined  first 
and  third  pressure  distributions,  several  maxima  and  minima  existed  on  the  disk. 
The  induced  velocity  was  directed  upward  at  the  front  and  sides  of  the  disk  and 
downward  at  the  rear  but  another  area  of  up  flow  existed  just  to  the  rear  of 
the  center  of  the  disk.  The  induced  velocity  at  the  center  of  the  disk  was 
equal  to  zero  and  was  independent  of  the  angle  of  incidence  of  the  rotor.  It 
was  found  that  the  local  peaks  in  the  upwash  could  be  removed  by  fairly  small 
alterations  of  the  load  distribution. 

At  large  distance  behind  the  rotor  the  vertical  component  of  the  induced 
velocity  was  directed  downwards  near  the  middle  and  upwards  on  both  sides  and 
was  symmetrical  with  respect  to  the  plane  « = o.  There  was  a singularity  . 
at  the  wake  boundary. 

In  reference  9,  the  Fourier  series  representation  of  the  downwash 
distributions  for  the  two  previously  mentioned  pressure  distributions  were 
given  in  terms  of  the  azimuth  angle,  rotor  radius,  thrust  coefficient,  flight 
path  velocity,  and  rotor  angle  of  incidence.  Only  pressure  distributions 
having  lateral  symmetry  were  considered.  This  Fourier  representation  makes  the 
calculation  of  the  blade  motion  more  convenient.  The  Fourier  series  could 
also  be  used  to  determine  numerical  values  of  the  downwash  distribution  over 
the  disk,  except  at  points  near  the  circumference  at  zero  incidence. 

The  authors  stated  in  reference  7 that  the  assumption  that  the  induced 
velocities  were  small  compared  to  the  flight  path  velocity  "implies  a serious 
limitation  in  the  applicability  of  the  results."  To  define  the  limits  of 
applicability,  the  following  criterion  was  established. 
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Ho  indication  was  given,  however,  as  to  the  "order  of  smallness  required  to 
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fulfill  this  criterion.  For  present  day  helicopters,  the  value  of  this  co- 
efficient falls  in  the  range  of  about  l/5  to  l/l2  with  the  cruising  speeds 
falling  toward  the  high  end  of  the  range  and  the  maximum  speed  felling  toward 
the  low  end. 

In  conclusion,  the  authors  stated  that  a sound  mathematical  theory  had 
been  constructed  for  the  induced  velocity  field  of  a rotor  within  the 
limitations  of  the  assumptions.  They  further  stated  that  there  should  be 
no  fundamental  difficulty  in  extending  the  theory  to  unsyrametrical  load 
distributions  which  occur  in  practice. 

An  attempt  was  made  by  P.  Brotherhood  and  W.  Stewart,  in  reference  12, to 
compare  the  variation  in  induced  velocity  along  a fore  and  aft  chord  of  the  - ! 

rotor  disk  as  calculated  by  the  theory  with  the  results  of  a full  scale  flight 
test  of  an  R-4B.  The  elliptical  load  distribution  was  used  in  the  theoretical 
calculation  for  both  0 and  15  degrees  angles  of  incidence  of  the  rotor.  It  was 
found  that  in  general,  this  distribution  resulted  in  values  which  were  higher 
(both  positively  and  negatively)  than  the  measured  values,  particularly  over 
the  aft  portion  of  the  disk.  The  measurements  were  made  for  yUy.  's  of  O.I38, 

0.167,  and  0.188.  The  agreement  was  somewhat  better  at  the  higher  flight 

1 

i 

speeds  than  at  the  lower,  indicating  possibly  that  the  limitations  of  the 
analysis  were  being  exceeded. 

In  reference  13  by  Daughaday  and  Kline  the  experimental  data  indicated 
the  presence  of  some  large  periodic  forces  exciting  the  first,  second,  and 

l 

third  bending  modes  of  the  rotor  blades  at  multip2.es  of  from  three  to  ten  times 
rotor  RPM.  These  results  were  in  complete  contradiction  with  the  theory  which 

; <1 

showed  that  the  exciting  forces  above  two  per  revolution  should  be  negligible. 

Since  the  theory  had  assumed  a uniform  downwash,  it  was  decided  to  use  Mangier' s 
downwash  analysis  to  calculate  the  generalized  forces.  The  combined  first  and 
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third  pressure  distributions  were  used  in  the  calculation.  The  results  of  this 
analysis,  according  to  reference  13> indicated  that  the  variation  in  induced 
velocity  was  probably  the  primary  source  of  the  higher  order  excitation.  It  was 
further  stated  that  the  theory  was  not  sufficiently  refined  to  be  used  ex- 
clusively in  rotor  blade  design,  in  that  the  designer  must  still  exercise 
considerable  judgement  in  the  selection  of  generalized  forces. 

Fail  and  Eyre,  in  reference  14,  made  some  downwash  measurements  behind  a 

12-foot  helicopter  rotor  in  a wind  tunnel  under  simulated  forward  flight  ! 

j 

conditions.  The  measurements  were  made  in  a vertical  plane  1.5  rotor  radii  I 

behind  the  rotor  center  at  four  spanwise  locations.  In  comparing  the  theory  and 
experiment,  an  allowance  must  be  made  for  the  displacement  of  the  rotor  wake 
since  the  first  order  theory  excluded  such  a displacement.  In  computing  the 
theory,  reasoning  seemed  to  indicate  that  the  first  or  elliptic  blade  spanwise 
pressure  distribution  would  most  nearly  approximate  the  loading  on  the  advancing  f 

r 

) 

blade,  while  the  combined  first  and  third  terms  of  the  Legendre  function,  which 

\ 

gave  a pressure  distribution  that  was  zero  at  the  rotor  center  and  tip,  would 
more  nearly  approximate  the  loading  on  the  retreating  blade.  In  general,  it 
was  found  that  the  experimental  data  approximately  corresponded  to  the  above 

it 

reasoning.  In  some  instances  the  agreements  between  theory  and  experiment  were 

f 

quite  good  if  allowance  was  made  for  the  wake  displacement.  In  other  cases  the  ! 

agreement  was  not  good.  There  appeared  to  be  no  correlation  between  jjy  ) aCy  , 
and  spanwise  location  and  the  agreement  or  disagreement  of  the  comparison.  In 
fact,  some  of  the  better  agreements  corresponded  to  flight  conditions  which 

violated  the  limitation  that  Cr  < < 1.  No  measurements  or  comparisbns  $ 

■i 

were  made  for  the  plane  of  the  rotor  blades. 
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Blade  Element  Theory 

Hie  simple  blade  element  theory  was  originally  Introduced  by  W.  Froude  in 
1878  (reference  15).  This  theory  considered  in  a crude  form,  the  forces 
experienced  by  the  blades  as  they  moved  through  the  air.  S.  Drzewiecki, 
beginning  about  1892  (references  16  & 17),  took  up  the  development  of  this 
concept  and  later  published  a complete  analysis  of  the  theory  (reference  18). 

The  analysis  was  based  on  the  premise  that  a propeller  blade  may  be  considered 
as  made  up  of  many  airfoil  elements  of  small  span  advancing  along  a helical 
path  determined  by  the  axial  velocity  and  the  rotational  velocity  of  the 
propeller.  Its  great  advantage  was  that  it  yielded  the  blade  geometry  but  had 
the  disadvantage  that  it  required  the  experimental  determination  of  the  airfoil 
characteristics  and  arbitrary  assumptions  as  to  the  effective  aspect  ratio.  At 
this  time  and  for  a considerable  period  thereafter  there  was  much  uncertainty 
as  to  what  airfoil  characteristics  should  be  used.  Drzewiecki  proposed  to 
determine  the  airfoil  characteristics  and  interference  effects  by  a series  of 

O 

tests  on  a number  of  special  propellers. 

F.  W.  Lanchester,  along  with  his  airfoil  theory,  developed  a blade  element 
theory  for  propellers  in  1907  (reference  19 ).  He  attempted  to  account  for  the 
mutual  interference  effects  by  an  anology  with  an  infinite  staggered  cascade 
of  airfoils. 

The  basic  assumption  of  these  theories  has  been  that  the  blade  element 
between  a radius  r and  r + A r of  a propeller  advancing  along  its  axis  with  a 
uniform  velocity  and  rotating  about  this  axis  with  a uniform  angular  velocity, 
was  considered  to  be  an  airfoil  element  of  the  same  cross-section  advancing 
through  the  air  with  a uniform  linear  velocity  which  was  the  resultant  of  the 
axial  and  rotational  velocities  and  at  an  angle  of  incidence,  cC  y . 
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Consequently,  the  aerodynamic  force  acting  on  the  "blade  element  was  the  same  as 
that  acting  on  the  aforementioned  fictitious  airfoil.  Then  the  total  force 
acting  on  the  blade  was  the  sum  of  that  acting  on  each  individual  blade  element. 
It  was  assumed  that  there  was  no  mutual  Interference  between  blade  elements 
except  as  the  interference  modifies  the  characteristics  of  the  blade  elements. 
According  to  these  primitive  blade  element  theories,  the  equations  representing 
the  complete  solution  of  the  behavior  of  a propeller  are 


4^  - TbcP  cv  cts  $ 

be**  f-  S>»-  <p  * ' Cjo  cos,  <pj  . 


where 
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Ci  = a(e-  <p) 
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These  equations  must  be  graphically  integrated  for  the  general  case. 
The  formula  for  the  efficiency  of  the  blade  element  is 

v,  _ / - e t*v  <p 


where  £ * 


‘ / + € cot  (p 

Solving  for  the  maximum  efficiency,  disclosed  the  fact 


that  according  to  the  blade  element  theory,  the  efficiency  of  a propeller 
decreased  as  the  diameter  increased  above  a certain  optimum  value  which  was  an 
absolute  contradiction  of  the  conclusion  reached  by  the  ideal  momentum  theory. 

Another  difficulty  was  that  these  theories  predicted  that  the  thrust  and 
torque  of  a given  propeller  at  a given  advance  ratio  should  vary  directly  as 
the  number  of  blades,  a fact  which  experimental  investigations  had  definitely 
disproved.  The  conclusion  which  these  discrepancies  indicated  was  that  the 
primitive  blade  element  theory  did  not  give  a complete  and  satisfactory  answer 
as  regards  to  the  behavior  of  the  propeller. 
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Combined  Momentum  and  Blade  Element  Theories 

It  was  pointed  out  in  the  previous  section  that  while  the  blade  element 
theory  predicted  the  observed  performance  of  propellers  in  a general  way,  it 
failed  to  give  accurate  numerical  results.  In  an  attempt  to  reconcile  the 
theory  with  experimental  results,  each  of  three  of  the  basic  assumptions  were 
critically  examined.  (See  reference  1,  page  215). 

The  first  assumption  considered  was  the  independence  of  blade  elements, 
i.e.  that  the  force  on  an  airfoil  element  was  not  affected  by  the  forces  on 
adjacent  elements  of  the  same  blade.  This  fact  could  not  be  established 
rigorously  but  the  effect  was.  believed  to  be  quite  small . However  in  order  to 
dispel  all  doubt,  a series  of  experiments  were  undertaken  by  C.  N.  H.  Lock, 
around  1924  (reference  20),  by  which  the  independence  of  blade  elements  was 
verified. 

The  next  assumption  to  be  considered  was  that  the  effective  velocity  of 
the  blade  through  the  air  was  the  resultant  of  the  propeller  axial  and  rotational 
velocity.  In  both  England  and  Germany,  it  was  proposed  that  the  increase  in 
axial  velocity  of  the  air  at  the  propeller  disk  as  predicted  by  the  momentum 
theory  should  be  included  in  the  blade  element  theory.  H..  Reissner,  in  1910, 
(reference  21)  was  the  first  who  tried  to  combine  the  momentum  theory  and  the 
blade  element  theory.  A.  Betz  in  1915;  (reference  22)  and  G.  de  Bothezat  about 
1918  (reference  23)  included  in  their  blade  element  theories,  this  increased 
axial  velocity  of  the  air  as  determined  by  the  ideal  momentum  theory.  In 
England,  however,  A.  Fage  and  H.  E,  Collins,  in  1917,  (reference  24)  proposed 

G 

that  an  empirical  estimate  of  the  increased  axial  velocity  of  the  air  be 
included.  This  latter  theory  was  extended  to  include  the  rotational  component 
by  L.  Bairstow  (reference  25).  Another  approach  to  the  problem  was  made  by 
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R.  McK.  Wood  and  H.  Glauert,  in  1918,  (reference  26).  They  used  Lanchester’s 
concept  (reference  19)  which  proposed  that  the  mutual  interference  of  propeller 
blades  was  analogous  to  that  of  a staggered  cascade  of  airfoils . Drzewiecki 
(reference  18)  maintained  that  the  increase  in  axial  velocity  of  the  air  at  the 
propeller  disk  as  determined  from  momentum  theory  was  only  an  average  of  a 
periodic  flow  and  therefore  should  not  be  used  in  estimating  the  force  experienced 
by  the  blade  element. 

The  third  assumption  was  concerned  with  what  airfoil  characteristics  should 
be  used.  All  of  the  previously  mentioned  blade  element  theories  suffered  from 
this  assumption.  Betz, in  reference  22 3 indicated  that  the  aspect  ratio  of  the 
propeller  blades  tended  to  be  infinite  but  depended  also  on  the  blade  shape; 

Fage  and  Collins,  in  reference  24,  used  an  aspect  ratio  of  6 and  attempted  to 
correct  the  theory  by  an  empirical  inflow  velocity.  Drzewiecki  and  de  Bothezat 
proposed  special  propeller  tests  to  determine  the  required  airfoil  characteristics. 
The  main  difficulty  was  that  the  variation  of  the  airfoil  characteristics  with 
aspect  ratio  was  still  undetermined  at  that  time,  both  from  a theoretical  and 
experimental  viewpoint. 

H.  Glauert  ir  1922,  (reference  27)  developed  a combined  axial  momentum  and 
blade  element  theory  in  which  he  indicated  that  to  be  consistent  with  the 
assumption  of  the  independence  of  blade  elements,  the  airfoil  characteristics 
must  be  taken  from  tests  in  which  the  airfoil  elements  experienced  no  interference 
from  adjacent  elements,  that  is,  infinite  aspect  ratio  tests.  This  analysis  was 
later  extended  to  the  case  of  the  helicopter  by  Glauert  (reference  28,  29,  and 
30),  C.  N.  H.  Lock,  (reference  31),  H.  B.  Squire  (reference  32)  and  many  others. 
These  references  gave  essentially  the  same  analysis  which  is  used  today  for 
helicopter  performance  and  design. 
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After  the  first  successful  flight  of  a helicopter,  it  became,  obvious  that 
these  theories  and  their  assumptions  were  not  adequate  for  predicting  the  various 
vibrations  and  blade  motions  encountered.  In  attempting  to  refine  the  theories, 
the  attack  was  directed  toward  finding  a theoretical  method  of  predicting  the 
variation  of  the  induced  velocity  across  the  rotor  disk  which  was  more  real- 
istic . 

Several  authors  have  made  use  of  the  combined  momentum  and  blade  element 
theory  in  this  field  with  some  success.  R.  S.  Ross  (reference  33)  set  up  the 
general  equations  determining  the  airflow  beneath  helicopter  rotors  using  these 
classical  concepts  as  applied  to  an  elementary  disk  annulus.  An  expression  was 
obtained  for  the  blade  span  loading  in  terms  of  the  local  induced  velocity  and 
the  local  lift-drag  ratio.  For  hovering  the  expression  became 
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In  order  to  solve  this  equation,  three  charts  were  prepared  as  follows:  1.  the 
particular  airfoil  characteristics,  and  £ versus  angle  of  attack;  2.  a 

plot  of  the  design  coefficient  - or  span  loading  versus  for  various 

lift-drag  ratios,  3.  a plot  of  Inflow  angle  f versus  velocity  ratio  for 

various  lift-drag  ratios.  Then  by  assuming  several  values  of  ci  , a value  of 
the  local  induced  velocity  could  be  determined  from  the  intersection  of  two  c, 
versus  • curves  for  the  particular  radius.  The  rotational  component,  u', 

could  then  be  determined  from  the  equation 
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Actual  velocity  measurements  were  made  with  a hot  wire  anemometer  underneath  a 
model  rotor  in  hovering  flight  and  compared  with  theory.  In  general  good  results 
were  obtained  over  the  inboard  center  portion  of  the  blade  for  both  velocity 
components  but  comparison  was  poor  toward  the  tips  and  close  to  the  root. 

The  design  formula  was  modified  to  account  for  horizontal  flight  with  the 
rotor  horizontal  and  then  further  modified  to  include  either  hovering,  vertical, 
or  horizontal  flight.  Although  velocity  measurements  were  also  made  for  the 
horizontal  flight  conditions,  there  was  no  comparison  with  theory.  Flight 
records  indicated  the  pulsations  or  periodicity  of  the  flow  that  occur  in 
horizontal  flight. 

W.  Castles  and  A.  L.  Ducoffe  (reference  34)  derived  an  expression  for  the 
distribution  of  the  induced  velocity  across  a rotor  disk  in  hovering  flight  as 
a function  of  the  thrust  at  that  radius  and  of  a semi-empi:  Leal  constant,  k, 
to  account  for  the  viscous  shearing  force  in  the  flow.  The  thrust  on  an 
elemental  annulus  as  determined  from  momentum  considerations  was  equated  to 
the  lift  on  the  portion  of  the  blades  lying  within  the  annulus  as  determined  by 
blade  element  considerations.  This  yielded  an  induced  velocity  distribution. 

V - 21  / Or-  k 

V "COS  ®Cj  \j  8 

where  «Cj,  is  the  wake  angle  of  each  annulus  stream  tube  at  the  disk.  Using  this 
value  of  the  induced  velocity,  the  induced  angle  of  attack  of  the  blade  was 
determined.  Then  using  the  relation  C- , = a(  — <py  ) and  the  determined 

relationship  <pr  = J ~z^£~  > the  local  lift  coefficient  for  each  blade 
element  could  be  determined.  The  thrust  and  torque  integral  equations  were  then 
set  up  for  any  twist,  taper,  or  planform.  These  integrated  equations  gave 
results  which  were  in  good  agreement  with  the  available  experimental  data.  The 
factor  k was  determined  empirically  from  experimental  data.  A similar  analyses 


22 


22 


was  made  for  hovering  within  ground  effect. 

As  pointed  out  in  the  previous  section  there  were  several  inconsistences 
and  contradiction  between  the  ideal  momentum  theories  and  the  blade  element 
theories.  One  inconsistency  was  that  in  direct  opposition  to  the  ideal 
momentum  theory  which  indicated  an  upper  limit  to  the  efficiency  of  a propeller 
depending  only  on  the  disk  loading  for  a given  advance  ratio,  the  blade  element 
theory  indicated  that  the  efficiency  approached  100$  as  the  drag  of  the  blade 
approached  zero.  Furthermore,  the  primitive  blade  element  theory  indicated 
that  the  efficiency  of  a propeller  decreased  as  the  diameter  increased  above  a 
certain  optimum  value  which  was  an  absolute  contradiction  of  the  ideal  momentum 
theory  which  indicated  that  the  efficiency  should  increase  with  diameter  for  a 
given  thrust.  The  rigorous  reconciliation  of  these  two  theories  awaited  the 
development  of  the  vortex  theory. 
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Vortex  Theory 


I 

l- 


!' 

i 

4 


ll 

t 

5* 


The  fundamental  bases  of  the  vortex  theory  were  the  Lanchester-Prandtl- 
Joukovskl  conception  of  the  vortex  system,  the  Kutta-Joukowski  theorem,  and 
Pran&tl's  airfoil  theory.  The  Kutta-Joukowski  theorem  stated  that  the  lift  of 
an  airfoil  was  determined  by  the  circulation  existing  about  its  contour. 

' . I 

That  is,  that  the  lift  per  unit  length  of  span  in  a two  dimensional  flow  was  j 

given  by  the  product  of  the  circulation  f about  its  contour,  the  density  p of 

the  fluid,  and  the  uniform  velocity  1/  of  the  airfoil  through  the  fluid,  namely  j 

I 
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The  concept  of  this  circulation  about  the  blade  led  to  the  concept  of  the  ! 

i 

1 

bound  vortex  or  lifting  line  and  as  a consequence  to  the  free  vortices  which 

are  shed  by  the  blade  and  pass  down  stream  in  a helical  pattern  forming  the 

slipstream.  Riis  type  of  slip-stream  configuration  was  apparently  first 

recognized  by  F.  W.  Lanchester  around  1900  and  was  later  published  in  1907  in 

reference  19*  N.  E.  Joukowski  investigated  the  induced  velocity  resulting 

from  this  type  of  vortex  system  in  1912  (references  5 and  35)  but  in  order  to  j 

reach  a solution  was  forced  to  assume  that  the  blades  were  lightly  loaded  and 

of  infinite  number.  This  resulted  in  a theory  identical  to  the  momentum  theory 

as  to  the  axial  and  rotational  induced  -velocities,  but  Joukowski  went  further 

and  suggested  using  an  infinite  aspect  ratio  cascade  of  airfoils  to  obtain 

the  required  airfoil  characteristics. 

A.  Betz,  in  1919  (reference  36), applied  Prandtl’s  elliptic  loading  airfoil 
theory  and  showed  that  a propeller  of  given  thrust  and  power  had  a minimum  loss 
of  energy  and  highest  efficiency  if  the  shed  vortex  sheets  after  an  initial 
limited  distortion,  moved  backward  as  a rigid  screw  surface  (i.e.  constant 
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axial  velocity).  The  application  of  this  theory  to  a propeller  was  not  simple, 
but  the  problem  was  capable  of  being  solved  if  the  propeller  was  assumed  to 
have  an  infinite  number  of  lightly  loaded,  frictionless  blades.  L.  Prandtl,  in 
an  addendum  to  reference  36,  indicated  an  approximate  method  for  correcting 
the  results  for  a propeller  with  a finite  number  of  blades.  In  this  case 
instead  of  a uniform  distribution  of  vortex  sheets  implying  no  gap,  there  was 
now  a system  of  vortex  sheets  with  a gap,  s , depending  upon  the  angle  <pt  of 
the  screw  surface  at  the  boundary  of  the  slipstream. 
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where  tan  <0  = and  b is  the  number  of  blades.  Whereas  before  there  was 

T*  -QJ? 

zero  slipstream  contraction  and  the  radial  velocity  component  was  negligible, 
there  is  now  an  appreciable  radial  velocity  at  the  boundary  of  the  slipstream 
as  the  air  attempts  to  flow  around  the  edges  of  the  vortex  sheets.  Prandtl 
replaced  these  vortex  sheets  by  a system  of  semi-infinite  plane  surfaces  of 
zero  thickness  and  gap,  S , equal  to  that  of  the  actual  vortex  sheet.  The  flow 
about  this  system  of  planes  was  determined  for  the  case  of  a uniform  motion 
downward  at  right  angles  to  the  planes.  Using  the  solution  to  this  system,  the 
flow  about  the  vortex  sheets  could  be  approximately  estimated.  A graph  of  the 
approximate  correction  factor  as  a function  of  the  number  of  blades  and  advance 
ratio  is  given  on  page  263,  reference  1.  S.  Goldstein  has  developed  a more 
accurate  analysis  of  the  problem  of  a finite  number  of  blades  by  the  use  of 
Bessel  functions  (reference  37).  This  was  a rigorous  solution  of  Beta's 
optimum  condition  but  assumed  zero  slipstream  contraction.  Hence  it  is  directly 
applicable  only  to  lightly  loaded  propeller.  Goldstein's  method  has  been 
extended  to  the  general  case  by  Lock,  reference  38.  These  theories  are 
applicable  to  a helicopter  rotor  only  for  the  special  cases  of  hovering  flight 
and  vertical  ascent.  A brief  summary  of  Lock's  theory  will  be  given  here. 
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A more  complete  form  of  the  derivation  may  be  found  in  reference  39. 

Goldstein  shoved  that  the  circulation  distribution  along  the  blades  of  a 
propeller  was  (according  to  his  assumptions) 

p _ 2 w' K »*  <f, 

1 ■ la 

where  2 W'  is  the  velocity  of  axial  translation  of  the  rigid  helicoidal 
surface;  K is  a non-dimensional  function  <p,  , the  ratio  , 

and  the  number  of  blades  b,  and  is  the  angle  that  the  helical  surfaces 

make  with  the  propeller  disk  plane.  It  was  then  shown  that 
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and  hence 


Lock  replaced  the  factor 
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represented  the  circulation  loss  along  the  blade.  Thus 
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by  a new  coefficient  X which 


Using  the  Kutta-Joukowski  theorem  that 


r * i^cV 


and  solving  for  the  axial  induced  velocity,  v , gave 
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The  rotational  component  was  then 
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The  expression  for  the  distribution  of  thrust  and  torque  along  a single  blade 


-M.  2 ± e cvl(cx  s/v  % + cdt  coi  ft) 


were 
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The  solution  of  these  equations  required  that  a value  of  <pf  he  assumed  for  each 
blade  station  considered  and  by  trial  and  error,  finding  the  distribution  of  <p, 
which  gives  the  desired  integrated  value  for  the  thrust.  This  is  at  best  a 
rather  complicated  solution  though  the  results  compare  very  well  with  experi- 
ments . 

M.  Knight  and  R.  Hefner, in  1937  (reference  40),  using  a vortex  system 
similar  to  that  of  Joukovski  (reference  5)  and  Glauert  (reference  27)  applied 
this  concept  to  the  lifting  rotor.  The  basic  assumptions  employed  were  the  same 
as  those  of  the  latter  two  references  and  were  as  follows:  the  number  of  blades 
was  taken  as  infinite;  the  induced  angles  were  small  so  that  the  sine  and 
tangent  were  equal  to  the  angle  itself  in  radians; and  the  rotational  and  radial 
components  of  velocity,  tip  losses,  and  slipstream  contraction  were  neglected. 
Each  blade  was  replaced  by  a rotating  lifting  line.  The  trailing  vortices  which 
sprang  from  the  tips  formed  a helix  in  space  while  the  vortices  from  the  root 
became  concentric  with  the  axis  of  rotation  and  were  neglected.  The  problem 
thus  was  reduced  to  finding  the  normal  component  of  the  velocity  which  was 
induced  at  the  rotor  disk  by  a cylindrical  surface  of  vorticity. 

In  order  to  simplify  the  analysis,  the  helical  pattern  of  vorticity  was 
broken  up  into  two  simpler  patterns.  One  was  composed  of  circular  vortex  rings 
and  the  other  of  axial  vortex  lines,  both  together  forming  a right  circular 
cylinder  of  vorticity  extending  from  the  rotor  disk  downward  to  infinity. 

However  in  this  analysis,  the  latter  pattern  was  neglected.  The  analysis  makes 
use  of  the  fact  that  the  potential  $ at  a point  P due  to  a closed  vortex 
ring  is  directly  proportional  to  the  product  of  the  circulation  P and  the 
solid  angle  subtended  by  the  vortex  ring  at  P . 
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The  solution  of  those  equations  required  that  a value  of  he  assumed  for  each 
blade  station  considered  and  by  trial  and  error,  finding  the  distribution  of  <p, 
which  (Ires  the  desired  integrated  value  for  the  thrust.  This  is  at  best  a 
rather  complicated  solution  though  the  results  compare  very  veil  with  experi- 
ments* 

M.  Knight  and  R.  Hefner,  in  1937  (reference  40),  using  a vortex  system 
similar  to  that  of  Joukovskl  (reference  5)  and  Glauert  (reference  27)  applied 

this  concept  to  the  lifting  rotor.  The  basic  assumptions  employed  were  the  same 
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as  those  of  the  latter  two  references  and  were  as  follows:  -the  number  of  blades 
was  taken  as  infinite;  the  induced  angles  were  small  so  that  the  sine  and 
tangent  were  equal  to  the  angle  itself  in  radians; and  the  rotational  and  radial 
components  of  velocity,-  tip  losses,  and  slipstream  contraction  vere  neglected, 
laoh  blade  was  replaced  by  a rotating  lifting  line.  The  trailing  vortices  which 
sprang  from  the  tips  formed  a helix  in  space  while  the  vortices  from  the  root 
became  concentric  with  the  axis  of  rotation  and  vere  neglected.  The  problem 
thus  was  reduced  to  finding  the  normal  component  of  the  velocity  which  was 
Induced  at  the  rotor  disk  by  a cylindrical  surface  of  vorticity. 

In  order  to  sinplify  the  analysis,  the  helical  pattern  of  vorticity  was 
broken  up  into  two  simpler  patterns.  One  was  composed  of  circular  vortex  rings 
and  the  other  of  axial  vortex  lines,  both  together  forming  a right  circular 
cylinder  of  vorticity  extending  from  the  rotor  disk  downward  to  infinity. 

However  in  this  analysis,  the  latter  pattern  was  neglected.  The  analysis  makes 
use  of  the  fact  that  the  potential  at.  a point  P due  to  a closed  vortex 

ring  is  directly  proportional  to  the  product  of  the  circulation  T ‘ and  the 
solid  angle  (0  subtended  by  the  vortex  ring  ”.t  P. 
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or  using  the  differential  equation  for  elementary  vortex  rings 

• 

This  equation  was  then  integrated  along  the  cylinder  from  i « 0 to  06  and 
then  differentiated  for  the  induced  velocity  at  the  disk.  It  was  then  shown 
that  for  a constant  circulation 


where 
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the  circulation  per  unit  length  along  the  cylinder.  It  was 


also  shown  that  the  induced  velocity  in  the  plane  of,  hut  outside  of  the  rotor 
disk  was  zero.  Furthermore,  it  was  shown  that  the  velocity  in  the  ultimate 
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thus  corroborating  the  momentum  result.  By  neglecting  higher  order  infinites- 
imals and  assuming  a variation  in  circulation  along  the  blade,  it  was  shown 
that  within  the  limitations  of  the  assumptions,  the  Independence  of  blade 
elements  held. 

These  results  were  used  to  derive  an  expression  for  the  inflow  angle  and 
hence  to  predict  the  lift  and  torque  on  a blade  by  means  of  the  usual  methods. 
The  resulting  thrust  and  torque  equations  compared  favorably  with  experiment. 

Knight  and  Hefner  also  made  a somewhat  similar  analysis  of  a lifting  air- 
screw in  ground  effect  (reference  4l).  All  the  previous  assumptions  were 
assumed  to  hold  except  the  Independence  of  blade  elements.  An  additional 
assumption  was  made  in  order  to  reach  an  approximate  solution.  This  was  that 
the  circulation  along  the  airscrew  blades  was  constant;  i.e.,  it  was  Independent 
of  both  the  blade  radius  and  the  distance  above  the  ground  plane.  The  vortex 
system  was  the  same  as  in  their  previous  paper,  except  that  the  vortex  cylinder 
extended  only  to  the  ground  plane.  The  effect  of  the  ground  plane  was 
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represented  by  placing  a second  cylindrical  vortex  sheet  of  equal  length  and 
strength  but  of  opposite  direction  at  the  end  of  the  first  cylinder.  Thus  the 
second  became  a mirror  image  of  the  first. 

The  authors  solved  for  the  value  of  the  induced  velocity  at  the  center  of 
the  rotor  and  at  the  tips  and  showed  that  the  solution  at  other  points  was 
extremely  difficult  and  tedious  due  to  the  presence  of  elliptic  integrals  of  the 
third  kind  for  which  no  tables  were  available.  Methods  of  approximation  were 
tried  and  Simpson’s  rule  yielded  satisfactory  results  for  values  of  ^ < 2, 
except  near  the  tips.  For  larger  values  of  or  as  X — » R,  this  method 

required  the  computation  of  a prohibitive  number  of  ordinates.  It  was  then 
shown  that  the  calculated  results  compared  satisfactorily  to  the  experimental 
data. 


R.  P.  Coleman,  A.  M.  Fe ingold,  and  C.  W.  Stempin  utilized  the  vortex 
theory  in  their  analysis  (reference  42)  of  the  longitudinal  variation  of  the 
induced  velocity  across  a rotor  disk  in  forward  flight.  The  study  was  under- 
taken in  order  to  find  a reasonable  explanation  for  the  vibration  encountered 
by  helicopters  in  slow  speed  forward  flight.  The  simplified  vortex  system  was 
also  used  as  a conception  of  the  rotor  wake.  In  this  instance  the  wake  pattern 
was  assumed  to  form  an  elliptic  cylinder  skewed  with  respect  to  the  rotor  axis 
at  an  angle  depending  upon  the  forward  velocity  and  the  induced  velocities.  It 
was  also  assumed  that  the  rotor  consisted  of  an  infinite  number  of  lightly  loaded 
blades  with  constant  circulation  and  that  the  helical  vortex  system  could  be 
replaced  by  a series  of  circular  vortex  rings  whose  plane  remained  parallel  to 
the  rotor  disk,  and  a series  of  axial  vortex  lines. 

In  order  to  find  the  wake  skew  angle  , the  induced  velocity  components 
parallel  and  perpendicular  to  the  wake  axis  were  assumed  to  be  constant.  Then 
by  taking  the  line  integral  of  the  velocity  around  suitable  paths,  it  was  shown 
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that  the  ratio  of  the  induced  velocity  component  perpendicular  to  the  wake  axis 
in  the  direction  of  forward  flight  to  the  induced  velocity  component  parallel 
to  the  wake  axis  was  equal  to  the  tangent  of  half  the  wake  angle. 

By  applying  the  Biot-Savart  law,  it  was  shown  that  the  velocity 
distribution  normal  to  the  disk  induced  by  the  assumed  vortex  wake  was  given 
by  a double  integral  equation.  This  equation  was  solved  at  the  center  of  the 
disk  and  the  value  obtained  was  then  shown  to  be  the  average  across  the  fore- 
and-aft  diameter.  Furthermore,  it  was  shown  that  this  value  was  half  the 
component  of  the  ultimate  wake  Induced  velocity  which  was  perpendicular  to  the 
rotor  disk.  The  fore-and-aft  rate  of  change  of  the  induced  velocity  was  found 
by  differentiating  within  the  integral  sign  and  evaluating  the  resultant 
integrals  at  the  center  of  the  disk.  It  was  found  that 

yc  s \ ft*  -f- 

vhsre  V.  « rate  of  change  of  induced  velocity  at  the  center  of  the  disk,  Ya  « 
average  value  of  fore-and-aft  Induced  velocity  at  the  center  of  the  disk,  and 
X » the  wake  skew  angle. 

The  results  were  combined  with  Glauert's  theory  of  reference  28  and  an 
expression  vas  arrived  at  for  the  induced  velocity  at  the  center  of  the  disk 
in  terms  of  the  flight  velocity  and  required  thrust.  A comparison  vas  made 
with  some  experimental  data  and  the  authors  concluded  that  this  comparison 
indicated  that  the  most  significant  factors  had  been  taken  into  account. 

In  reference  12,  the  results  of  this  analysis  were  also  compared  with  the 
experimentally  determined  induced  velocities.  In  this  case  the  measured  values 
of  the  Induced  velocities  at  the  center  of  the  disk  were  somewhat  lower  than 
the  calculated  values  for  all  three  flight  speeds.  The  experimentally  determined 
slopes  of  the  induced  velocity  variation  were  all  higher  than  that  yielded  by 
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the  analysis.  The  theory  appeared  to  yield  a better  comparison  at  the  lower 
flight  speeds.  No  definite  conclusions  were  reached,  however,  since  the 
experimental  induced  velocities  were  not  determined  for  the  fore-and-aft 
diameter,  and  since  some  lateral  asymmetry  was  supposed  to  exist. 

In  reference  13,  as  stated  in, a previous  section,  large  periodic  forces 

were  found  to  exist  on  rotor  blades  in  forward  flight  which  excited  the  first, 

second,  and  third  bending  modes  of  the  rotor  blades  at  frequencies  of  from 

three  to  ten  times  rotor  frequency.  As  before  an  effort  was  made  to 

theoretically  predict  these  forces,  and  the  linear  downwash  distribution  of 

Coleman,  Feingold,  and  Stempin  was  incorporated  in  the  theory  for  finding  the 

generalised' fprces . It  was  found  thatt  the  goading  terms  which  were  directly 

affected  by  this  type  variation  were  the:  first  harmonic  cosine  loading  and  the 

second  harmonic  sine  loading.  The  remaining  higher  order  harmonic  loadings 

were  proportional  to  blade  flapping  coefficients.  It  was  shown  that  this 

linear  induced  velocity  distribution  had  a minor  effect  on  all  blade  flapping 

coefficients  except  the  first  harmonic  sine  coefficient.  A comparison  was  made 

of  the  theoretical  generalized  forces  and  the  test  data.  It  was  found  that  fair 
» 1’  • ■ ' 
agreement  was  obtained  at  high  tip  speed  ratios,  but  the  analysis  was  not 

capable  of  explaining  the  peaks  in  the  generalized  forces  in  the  transition 

region. 

Drees  (reference  43),  by  using  the  vortex  concept  of  the  lifting  rotor, 
has  developed  a theory  which  is  applicable  to  all  flight  conditions  but  was 
still  limited  to  a great  extent  by  the  necessary  assumptions  for  the  solution 
of  the  problem.  In  his  paper,  the  usual  assumptions  were  made:  i.e. 
contraction  and  rotation  of  slipstream  was  neglected,  number  cf  rotor  blades 
was  infinite,  circulation  was  constant  along  the  blade  radius,  and  the  angle 
between  the  slipstream  and  the  rotor  disk  at  each  point  of  the  disk  was  constant. 
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Although  It  vas  assumed  that  the  circulation  along  the  blade  was  constant,  it 
vas  quite  obvious  that  in  forward  flight,  the  circulation  must  vary  with  the 
azimuth  angle  . It  vas  then  assumed  that  the  circulation  could  be  written 

as 

r = r0  - r,  **  + 

Values  of  To  and  P,  , were  found  by  equating  the  constant  portion  of  the 
lift  to  the  weight  and  by  setting  the  condition  that  the  lift  moment  about  the 
flapping  hinge  must  be  constant  and  independent  of  ^ . Then  by  neglecting 
second  and  higher  harmonics  it  was  shown  that 

r° = 

r,  - 

The  axial  velocity  was  taken  as  being  composed  of  the  following  components : 
the  component  of  the  flight  velocity  perpendicular  to  the  tip  path  plane ; the 
component  due  to  the  bound  vortices  around  the  blades;  the  component  from  the 
free  vortices  moving  with  the  slipstream  and  springing  from  the  blade  tips 
only,  since  P is  assumed  constant  along  the  radius;  and  the  component  from  the 
free  radial  vortices  influencing  the  axial  flow  which  compensates  for  the 
variation  of  T with  respect  to  Ip  . The  induced  velocities  were  calculated 
only  for  the  X and  Y axes  at  the  3/U  radius  station  and  at  the  rotor  center. 

A more  accurate  solution  was  considered  impractical  in  use  and  not  in  accordance 
with  the  other  simplifications.  For  performance  calculations,  a certain 
velocity  distribution  vas  assumed  to  exist  through  these  points  as  a function  of 
the  radius  and/or  as  a first  harmonic  function  of  the  azimuth  angle,  & . 

Terms  of  higher  than  the  first  order  were  neglected. 

The  induced  velocity  of  the  bound  vortices  was  arrived  at  by  means  of  the 
following  reasoning  and  by  the  usual  mathematical  solution  for  a finite  vortex 
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line.  It  was  reasoned  that,  since  /"*  was  symmetrical  to  the  Y axis,  it  could 
have  no  influence  on  the  axial  velocity  across  a lateral  diameter  of  the  disk 
and,  for  the  same  reason,  could  have  no  influence  on  the  axial  velocity 

along  the  X axis  or  longitudinal  diameter  of  the  rotor  disk.  Thus  only  /") 
gave  an  axial  induced  velocity  component  along  the X axis.  The  solution  of 
the  mathematical  relationships  involved  complete  elliptic  integral  equations. 
Assuming  the  component  of  the  induced  velocities  due  to  the  hound  vortices  to 
he  a cosine  function,  these  equations  yielded  for  the  3/4  radius  station 
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The  component  of  the  induced  velocity  due  to  the  circulation 
component  of  the  free  tip  vortices  was  found  in  the  usual  manner  hy  assuming 
that  the  shed  vortices  formed  closed  vortex  rings  parallel  to  the  tip  path 
plane.  Two  integral  equations  were  arrived  at  for  the  velocities  along  the 
longitudinal  and  lateral  rotor  diameter,  hut  the  integration  while  possible, 
was  very  complicated.  These  equations  were,  however,  integrated  for  the  rotor 
center  and  yielded 
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where  is  the  angle  between  the  slipstream  and  the  perpendicular  to  the 

tip-path  plane.  When  the  equations  were  solved  for  oC4  =0  , they  yielded 
an  axial  induced  velocity  which  was  constant  for  each  point  of  the  rotor  disk. 
This  was  in  agreement  with  the  results  of  previous  papers,  which  showed  that 
the  assumption  of  constant  circulation  along  the  blades  gave  the  ideal  axial 
velocity  distribution  in  hovering.  The  two  integral  equations  were  also 
numerically  Integrated  for  the  3/4  radius  stations  on  both  longitudinal  and 
lateral  diameters.  Then  assuming  a cosine  variation,  the  result  was  given  as 
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or  in  another  manner 
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The  determination  of  the  axial  induced  velocity  due  to  the  free  radial- 
vortices  and  tip-vortices  was  more  difficult.  Since  P varied  cyclicly,  a 
free  axial  vortex  of  7 Jr  was  released  and  the  free  tip  vortex  also  varied 
t Jr  . These  vortices  formed  a surface.  Using  these  actual  distributions 
of  circulation  resulted  in  a very  complicated  calculation  for  the  component  at 
the  3 A radius  station.  In  order  to  arrive  at  a solution,  it  was  necessary  to 
resort  to  a simple  analogy.  The  radial  and  tip  vortices  were  taken  together 
and  assumed  to  form  two  circular  cylinders.  That  is,  the  actual  circulation 
distribution  was  replaced  by  the  vortex  system  of  two  hypothetical  rotors. 

These  rotors  had  their  axis  on  the  lateral  diameter  of  the  actual  disk.  The 
cylinder  at  ’jr  ■ 90°  had  a strength  of  — , while  the  cylinder  at 

'f'  * 270°  had  a strength  of  + f~j  . These  systems  condensate  on  the 
longitudinal  diameter  and  hence  cause  no  variation  in  the  Induced  velocity. 

Then  in  a maimer  similar  to  that  of  the  free  constant  tip -vortices,  the  induced 
velocity  component  at  the  3/4  radius  station  and  assuming  a sine  variation  was 
shown  to  be 
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All  the  above  components  were  added  together  to  yield  the  total  axial 
induced  velocity.  The  equations  were  then  non-dimensionalized  but  were  not 
valid  for  ideal  vertical  autorotation  when  the  airflow  through  the  rotor  disk 
approached  zero.  The  theory  was  then  changed  in  a semi -empirical  manner  to 

V 

Include  this  flight  condition.  A formula  was  arrived  at  which  was  valid  for 


all  flight  conditions  but  had  three  constant  coefficients  which  must  be 
determined  from  flight  and  wind-tunnel  tests.  A correction  for  the  deviation 
of  the  actual  from  the  ideal  uniform  induced  axial  velocity  was  included.  A 
chart  was  then  presented  which  is  valid  for  all  kinds  of  one-rotor  helicopters 
and  autogyros  and  all  normal  speed  combinations.  The  paper  concluded  with  the 
equations  for  the  determinations  of  «C(  , the  attitude  angle  of  the  fuselage 
and  SL,s  , the  angle  between  the  rotor  axis  and  the  perpendicular  to  the 
tip-path-plane.  Applications  of  the  theory  to  several  problems  were  given. 
Performance  calculations  were  made  for  the  Sikorsky  S-51  and  the  theory  pre- 
dicted the  experimental  data  to  a good  degree  of  accuracy. 

W.  Castles,  Jr.  and  J.  H.  De  Leeuv,  in  reference  44,  presented  a method 
for  computing  the  approximate  values  of  the  normal  component  of  the  induced 
velocity  at  points  in  the  flow  field  of  a lifting  rotor.  In  this  paper  the 
assumption  was  also  made  that  the  slipstream  of  the  rotor  could  be  considered 
as  a uniform  continuous  distribution  of  vortex  rings  of  infinitesimal,  strength. 
These  rings  were  assumed  to  lie  in  planes  parallel  to  the  tip  path  plane  and  to 
fora  a straight  elliptic  cylinder  extending  from  the  disk  to  infinity.  Then 
by  using  the  stream  function  at  a point  P in  the  flow  field  of  a vortex  ring, 
expressions  were  arrived  at  for  the  axial  and  radial  components  of  the  induced 
velocity.  These  expressions  were  rather  complicated  since  they  involved  the 
complete  elliptic  integrals  of  the  first  and  second  kinds  and  were  given  as 
follows : 
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where  the  terms  within  the  brackets  are  functions  of  the  complete  elliptic 
integrals  and  non-dimensional  distances  of  the  point  P from  certain  references. 
These  terms  are  defined  in  the  original  paper. 
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It  nay  be  seen  that  the  egressions  become  indeterminant  for  x « o (rotor 
center) . It  had  been  previously  shown  that  for  this  case . 
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The  above  expressions  were  rewritten  in  non-dimensional  form.  The  factor. 


the  non-dinensional  normal  component  of  the  induced  velocity  in  the 


vicinity  of  a vortex  ring  appeared  on  the  left  side  of  the  expression  and  the 
right  side  of  the  expression  became  a function  of  the  geometry  of  the  problem; 
i.e.  the  distances  of  the  point  P from  the  center  of  the  ring.  Numerical 
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values  of  this  factor 


were  calculated  and  tabulated  for  various  non- 


dimensional  axial  and  radial  distances  from  the  center  cf  the  ring.  The  intervals 
were  so  chosen  that  numerical  integration  by  Simpson's  rule  could  be  accomp- 
lished. Then  using  the  resulting  table,  the  normal  component  of  induced  velocity 
at  a point  P in  the  vicinity  of  a rotor  whose  wake  was  composed  of  these  rings 
could  be  found  by  adding  up  the  contributions  of  each  vortex  ring  covered  by 
the  table.  This  would  account  for  about  95$  of  the  total  induced  velocity 
component  at  the  center  of  the  rotor  and  for  a large  part  of  the  component  for 
most  points  considered  by  the  paper.  The  contribution  of  the  rings  outside  the 
scope  of  the  table  were  summed  by  using  an  approximate  integral  equation  capable 
of  solution  which  was  said  to  introduce  small  error.  Thus  an  expression  was 
arrived  at  for  the  contributions  of  the  vortex  rings  from  the  limits  of  the 
table  to  infinity.  This  expression  was  a function  of  the  geometrical  position 
of  the  point  under  consideration  and  the  wake  angle  of  the  slip-stream  at  the 
rotor. 

The  results  were  presented  in  the  form  of  tables  and  graphs  of  the  ratio 
of  the  normal  induced  velocity  component  at  point  P to  the  normal  component 
of  induced  velocity  at  the  center  of  the  rotor. 
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A method  was  then  given  for  the  determination  of  the  mean  value  of  the  normal 
component  of  induced  velocity  over  the  front  and  rear  rotors  of  a tandem 
helicopter;  for  the  determination  of  the  longitudinal  variation  of  the  normal 
component  of  the  induced  velocity  over  the  front  and  rear  rotors  of  a tandem 
helicopter  and  for  the  determination  of  the  induced  flow  angle  at  a horizontal 
tall  plane. 

In  conclusion,  the  author?  of  the  paper  stated  that  for  the  high  speed 
forward  flight  condition,  the  assumption  that  the  vortex  . ^ngs  remain  parallel 
to  the  tip  path  plane  was  the  only  one  likely  to  affect  the  engineering 
accuracy  of  the  results.  Whereas  for  the  slower  speed  case,  the  initial 
assumptions  as  to  the  vortex  ring  spacing  and  circulation  distribution  along 
the  blades  would  affect  the  accuracy  of  the  calculations.  They  therefore 
advised  exercising  caution  in  applying  the  theory  to  points  on  or  close  to 
the  disk  of  a specific  rotor  at  low  forward  flight  velocities.  No 
comparison  with  experiment  was  made . 
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Flov  Visualization  Studies 

There  have  been  many  experimental  investigations  which  have  studied  the 
air  flov  through  and  around  both  propellers  and  helicopters.  Most  of  these 
have  dealt  with  models,  though  a few  have  used  full  scale  helicopters  under 
actual  flight  conditions.  The  purpose  vas  to  gain  a more  accurate  knowledge 
of  the  flows,  with  the  thought  that  if  the  actual  flov  were  understood, 
perhaps  the  theory  could  be  improved  either  by  incorporating  simplifying 
assumptions  in  the  mathematics  or  by  devising  an  empirical  method.  The  results 
were  also  used  as  a check  on  the  applicability  of  existing  theories. 

There  have  been  four  methods  of  flow  visualization  employed:  cotton  tufts, 
smoke,  balsa  dust,  and  hot  wire  and  spark  shadowgraphs.  From  the  available 
experimental  results,  the  following  general  conclusions  may  be  drawn.  Cf  the  _ 
four,  the  cotton  tufts  seem  to  be  of  least  value,  since  they  suffer  a great 
deal  from  gravity  and  flutter,  particularly  in  the  rotor  wake.  Smoke  photo- 
graphs appear  to  yield  better  results  for  streamline  study,  while  balsa-dust 
photographs  appear  to  yield  a better  overall  picture  of  the  flow.  The  hot 
wire  and  spark  shadowgraphs  appear  to  be  the  better  method  for  studying  vortex 
and  vake  helix  patterns  because  the  heated  filaments  of  air  do  not  appear  to 
dissipate  too  quickly.  Some  examples  of  these  methods  as  applied  to 
propellers  and  rotors  are  reviewed  in  the  following  paragraphs. 

Lock  and  Tovnend,  in  reference  45,  investigated  the  flow  around  a model 
propeller  working  in  water  in  the  "vortex  ring  state."  The  model  propeller 
was  of  extremely  small  scale  having  a diameter  of  one  inch.  Ten  photographs 
were  presented  including  views  of  the  flow  conditions  from  the  normal  working 
state  through  the  static  thrust  condition  to  the  vortex  ring  and  windmill 
brake  state.  Lock,  in  reference  46,  repeated  these  tests  using  a three  foot 
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diameter  airscrew  in  a windtunnel  with  cotton  tufts  as  the  means  for  flow 
visualization.  Twelve  photographs  were  presented  covering  the  various  states 
as  listed  above.  The  photographs  clearly  showed  the  extreme  turbulence  which 
exists  about  an  airscrew  in  the  vortex  ring  state. 

Townend,  in  reference  47/  used  hot  wire  and  spark  shadowgraphs  in  a study 
of  the  periodic  airflow  through  a propeller.  The  use  of  the  spark  shadowgraphs 
enabled  the  velocity  and  the  direction  of  flow  to  be  mapped  out.  The 
shadowgraphs  of  the  hot  wire  tests  very  clearly  showed  the  radial  flow  that 
exists  in  the  slipstream.  The  heated  air  filaments  were  broken  by  the  blade 
and  this  break  passed  downstream  as  a discontinuity,  the  particle  that  passed 
just  under  the  blade  moved  outward  while  the  adjacent  particle  which  had 
passed  just  over  the  blade  moved  inward.  This  discontinuity  represented  the 
helicoidal  vortex  sheet  left  behind  by  the  blade.  The  filaments  also  very 
clearly  showed  the  tip  vortices  which  were  shed  and  their  spacing.  The  centers 
of  these  "eddies"  as  they  appeared  in  the  shadowgraphs  represented  the  slip- 
stream boundary.  The  heated  dots  formed  by  the  injected  sparks  were  intro- 
duced near  the  blade  tips.  The  resulting  shadowgraphs  showed  very  clearly 
the  formation  of  the  trailing  tip  vortices.  The  high  velocities  which  existed 
in  this  region  were  indicated  by  the  elongation  of  the  dot  shadows.  These 
photographs  are  an  excellent  though  brief  study  of  the  periodic  flow  existing 
behind  an  airscrew. 

Brotherhood,  in  references  48  and  49,  and  Brotherhood  and  Stewart,  in 
reference  12,  investigated  the  air-flow  through  a full  scale  helicopter  rotor 
in  various  free  flight  conditions.  Velocity  measurements  were  made  and  pictures 
were  taken  of  the  various  flow  patterns  formed  by  smoke  streamers.  Reference 
48  is  an  investigation  of  the  air-flow  through  a helicopter  rotor  when  hovering 
both  in  and  out  of  ground  effect.  The  smoke  visualization  of  the  streamlines 
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was  used  to  refer  back  to  their  appropriate  disk  positions  the  velocity  measure- 
ments  which  were  taken  in  the  fully  developed  slipstream.  The  experimental 
measurements  agreed  very  well  with  that  predicted  by  propeller  strip  theory  of 
reference  1,  if  tip  losses  were  taken  into  account.  The  smoke  filaments  gave 
a good  indication  of  the  streamline  pattern  in  the  slipstream.  Reference  49 
was  an  investigation  of  the  air-flow  through  a rotor  in  vertical  descent. 

Smoke  filaments  were  Introduced  below  the  rotor  at  several  radii.  The  pictures 
of  the  smoke  pattern  in  vertical  descent  clearly  indicated  the  various  types 
of  flow  associated  with  the  rotor  as  it  passed  through  the  vortex-ring  state 
into  the  windmill-brake  state  of  operation.  There  was  very  good  agreement 
between  the  rates  of  descent  corresponding  to  autorotation  as  obtained  from  the 
smoke  photographs  and  from  power  considerations.  Reference  12  has  been 
previously  discussed.  It  presented  smoke  photographs  of  the  flow  through  a 
rotor  in  forward  flight . The  helicopter  was  flown  behind  an  aircraft  which 
was  trailing  smoke  generators  suspended  from  a long  wire.  The  smoke  trails 
passed  through  the  rotor  and  were  photographed  from  the  side  by  another  air- 
craft. The  resulting  photographs  indicated  the  increase  in  induced  velocity 
from  the  front  to  the  rear  of  the  disk.  The  results  were  in  reasonable  agreement 
with  theoretical  predictions. 

Taylor,  in  reference  50;  developed  and  illustrated  a balsa-dust  technique 
for  the  visualization  of  air-flow  patterns  through  model  helicopter  rotors  both 
in  steady-state  and  transient-flow  conditions . Some  very  striking  photographs 
were  obtained,  particularly  of  the  starting  vortex  shed  from  the  blade  tips 
while  the  rotor  was  brought  up  to  speed  from  rest,  both  for  single  and  coaxial 
rotor  configurations.  Photographs  were  also  given  of  single,  coaxial,  and 
biaxial  (with  varying  degrees  of  overlap)  rotor  configurations  in  and  out  of 
ground  effect.  This  method  illustrated  the  overall  flow  picture  very  well. 
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Drees  and  Hendal,  in  reference  $1,  photographed  the  flow  of  smoke  filaments 
through  and  around  a model  helicopter  mounted  in  a wind  tunnel  under  various 
flight  conditions.  A method  was  devised  for  introducing  smoke  at  the  exact  local 
velocity.  The  smoke  could  be  turned  on  and  off  instantaneously,  was  non- 
corrosive,  and  was  easily  regulated.  Photographs  of  the  smoke  flow  were  made 
for  hovering  flight  and  various  rates  of  vertical  descent  through  the  vortex- 

O. 

ring  state  into  the  windmill -brake  state.  In  general  the  turbulence  was  so 
great  as  to  completely  dissipate  the  smoke  filaments  after  they  had  passed 
through  the  rotor.  An  interesting  series  of  photographs  showed  a marked 
periodicity  of  airflow  pattern  for  descending  flight  at  lew  forward  speed,  the 
helicopter  being  in  the  vortex-ring  state.  This  corresponds  to  the  region  of 
roughness.  The  model  helicopter  rotor  was  observed  to  "tumble"  regularly  with 
the  same  period  as  the  shedding  of  the  vortex  ring.  At  higher  forward  speeds, 
the  phenomenon  disappeared.  For  autorotation  in  forward  flight,  the  photo- 
graphs showed  that  there  was  up  flow  in  the  front  part  of  the  disk,  but  that 
near  the  rear,  there  was  a region  of  downflow. 

Meyer  and  Falabella,  in  reference  52,  have  measured  the  aerodynamic 
loading  on  a model  helicopter  blade  in  forward  flight  by  recording  the  pressure 
variations  by  means  of  pressure  taps,  positioned  both  spanwise  and  chordwise 
on  the  blade.  The  results  have  been  plotted  as  curves  of  constant  aerodynamic 
loading  versus  position  on  the  rotor  disk.  This  reference  would  yield 
excellent  data  for  checking  any  future  theory. 


4l 


Discussion 

As  may  "be  seen  from  the  reviews  presented  in  the  preceeding  sections, 
simplifications  of  the  physical  picture  were  necessary  in  order  to  obtain 
solutions  for  the  induced  velocity  field.  These  assumptions  and  their  effects 
will  be  discussed  below. 

The  primary  and  perhaps  the  most  universal  assumption  is  that  the  rotor  is 
made  up  of  an  infinite  number  of  blades.  In  the  vortex  theory  this  implies  zero 
slipstream  contraction  since  the  resulting  wake  boundary  may  then  be  taken  as  a 
cylindrical  sheet  of  vorticity  for  which  there  can  be  no  radial  flow  within  the 
wake.  In  these  cases  the  additional  assumptions  are  usually  made  that  the 
circulation  along  the  blades  is  constant  and  that  the  effect  of  the  vortices 
shed  off  the  blade  root  may  be  neglected  except  insofar  as  it  results  in  a 
rotational  component  of  induced  velocity  about  the  wake  axis.  In  the  momentum 
theory,  this  assumption  implies  that  the  rotor  may  be  replaced  by  a pressure 
step,  the  distribution  of  which  must  either  be  assumed  or  solved  for  with  the 
aid  of  other  simplifying  assumptions  while  the  velocity  remains  continuous 
through  the  disk.  An  extension  of  this  concept  is  that  of  the  actuator  disk 
in  which  the  rotor  is  replaced  by  a pair  of  Infinitely  bladed,  contra-rotating 
tandem  rotors  so  designed  that  the  rotational  component  of  induced  velocity  of 
the  front  rotor  is  exactly  cancelled  by  the  rear  rotor.  The  result  of  the 
assumption  of  an  infinite  number  of  blades  yields  an  induced  velocity  variation 
which  is  not  periodic.  For  fast-moving  multi-bladed  propellers  and  rotors  this 
would  appear  to  be  a good  assumption  but  for  large-diameter,  slow-moving  rotors 
with  only  one  or  several  blades,  the  assumption  would  appear  to  be  rather  poor 
even  for  performance  estimation. 

The  second  usual  assumption  is  that  the  blades  are  lightly  loaded.  This  may 
also  imply  several  things . One  is  that  the  slipstream  contraction  may  again  be 
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neglected  which  means  that  the  inflow  velocity  must  be  very  small  compared  to  the 
rotational  velocity  component  of  the  blade  section.  Another  is  that  in  forward 
flight,  the  velocity  along  the  flight  path  is  large  compared  to  the  induced 
velocities  and  hence  the  theory  may  be  linearized.  This  assumption  may  result 
in  a reasonable  approximation  if  care  is  taken  in  interpreting  the  limitation 
of  smallness.  At  present  this  limitation  is  not  clearly  defined. 

The  assunptions  that  the  slipstream  contraction  and  rotation  may  be 
neglected  while  in  error  physically  appear  to  have  little  or  no  effect  on  the 
results  of  the  vortex  theory  as  compared  to  performance  data.  However,  the 
former  is  obviously  incompatible  with  the  momentum  theory  while  the  latter 
probably  has  a small  effect. 

The  assumption  of  a constant  circulation  along  the  blade  span  axis  has 
been  a necessary  assumption  in  most  vortex  theories.  (The  exception  being 
Goldstein's  theory.)  This  distribution  is  theoretically  obtainable  by  the 
proper  twist  and/or  taper  of  the  blade  but  this  is  not  practicable.  In  practice 
the  actual  circulation  distribution  is  quite  different  from  the  assumed  constant 
value. 

It  has  been  shown  experimentally  by  Lock  in  reference  20  that  the 
assumption  of  the  independence  of  blade  elements  in  the  blade  element  theories 
is  a reasonable  approximation  and  yields  good  results. 

Goldstein's  theory  in  reference  37  and  as  modified  by  Lock  in  reference  38 
employed  only  one  of  the  assunptions  stated  above  which  was  that  the  blades  must 
be  lightly  loaded,  though  this  later  was  shown  to  be  an  unnecessary  condition. 

The  basic  assumption  of  the  theory  was  that  the  wake  was  pictured  as  composed 
of  a rigid  helicoidal  surface  (or  co-axial  surfaces)  of  infinite  length  but 
finite  radius  moving  with  a uniform  axial  velocity  which  was  small  compared  to 
the  rotor  blade  tip  speed.  Then  the  velocity  field  of  the  vortex  system  at  a 
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great  distance  behind  the  rotor  was  equivalent  to  the  potential  field  of  these 


rigid  helicoidal  surfaces.  The  conditions  existing  in  the  wake  were  then  related 


to  the  rotor  planform.  The  solution  is,  however,  only  rigorously  correct  for 


rotors  whose  distribution  of  circulation  along  the  blades  will  yield  a flow  in 


the  wake  that  is  identical  with  the  potential  flow  of  such  a set  of  rigid 


equidistant  co-axial  helicoidal  surfaces.  It  appears  that  the  theory  may  be 


applied  with  small  error  to  the  general  case  of  a hovering  or  vertically 


ascending  rotor.  This  is  the  only  theory,  with  the  exception  of  Prandtl’s 


tip  loss  correction  (reference  36),  that  takes  into  account  the  interference 


effects  of  a finite  number  of  blades. 


In  all  cases  the  working  medium  was  assumed  to  be  an  incompressible. 


inviscid  gas. 


In  the  main,  these  assumptions  which  have  been  discussed  appear  to  be 


reasonably  well  founded  because  the  resulting  performance  theory  predicts  the 


overall  experimental  data  to  a good  degree  of  accuracy.  However,  when  these 


theories  are  applied  to  the  calculation  of  the  dynamic  air  loadings  on  the 


blade  for  blade  motion,  vibration,  and  stress  analyses,  it  becomes  apparent 


that  they  are  woefully  inadequate  (for  instance,  the  efforts  of  Doughaday  and 
Kline  in  reference  13).  It  would  therefore  appear  to  be  worthwhile  to  attempt 


to  extend  the  present  theories  or  to  find  another  solution  to  the  problem 


which  is  not  as  greatly  limited  in  its  applicability.  One  of  the  primary 


objectives  of  such  an  attempt  would  be  the  incorporation  of  the  effect  of  the 


periodicity  of  flow  and  the  interference  effects  of  a small  number  of  blades. 


The  neglect  of  these  effects  may  be  the  cruse  of  the  major  discrepancies 


appearing  in  the  present  theories, 


The  search  for  a general  solution  has  been  almost  completely  thwarted  by 


the  lack  of  suitable  mathematical  tools.  This  deficiency  has  made  it  necessary 
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to  so  simplify  the  true  physical  picture  that  many  of  the  practical  uses  vhich 
could  have  been  derived  from  the  results  have  largely  been  lost.  It  is  obvious 
from  the  reviewed  papers  that  men  of  stature  have  grappled  with  the  problem  and 
have  met  only  limited  success,  not  through  lack  of  effort  or  intelligence,  but 
because  of  the  limits  of  mathematical  knowledge.  The  author  of  this  paper  does 
not  presume  to  have  the  ability  to  extend  these  mathematical  limits.  However, 
with  these  thoughts  in  mind,  a program  of  investigation  has  been  laid  out 
which  should  result  in  some  further  degree  of  success. 

It  would  seem  that  the  possibility  of  practically  extending  the  knowledge 
by  using  the  mathematical  artifices  and  physical  simplifications  of  the 
reviewed  papers  have  been  largely  exhausted.  Further  advances  along  these 
lines  would  appear  to  require  prohibitive  amounts  of  graphical  integrations  or 
numerical  calculations  though  electronic  computers  may  make  this  approach 
feasible. 

It  is  therefore  proposed  that  a search  for  a solution  be  instituted 
which  will  be  guided  along  two  lines  of  endeavor.  The  first  of  these  would 
be  a systematic  search  of  the  more  recent  papers  and  journals  in  the  fields 
of  mathematics  and  the  physical  sciences  for  a new  or  different  mathematical 
method  or  artifice.  The  second  approach  would  be  a flow  visualization  study 
of  the  velocity  field  about  a helicopter  rotor  in  the  various  flight  conditions. 
The  purposes  of  the  latter  study  would  be  to  indicate  what  simplifications  may 
be  made  vhich  would  not  too  greatly  limit  the  applicability  of  a solution  to 
practical  problems.  It  is  thought  that  these  two  methods  of  approach  would  offer 
some  hope  of  finding  a more  satisfactory  solution  to  the  problem  of  the  induced 
velocity  field  about  a helicopter  rotor. 
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Conclusions  8c  Recommendations 

The  rotor  induced  velocity  field  theories  that  have  been  reviewed  herein 
have  made  definite  and  worthwhile  contributions  to  the  understanding  of  the 
flow  through  helicopter  rotors.  However  they  have  in  the  main  been  either 
over-siaplified,  limited  by  their  assumptions  to  special  flight  conditions,  or 
are  rather  long  and  tedious  in  application.  In  general  the  problem  of  the 
periodicity  of  flow  has  been  neglected  entirely.  It  is  therefore  recommended 
that  a program  of  study  be  initiated  which  will  be  directed  toward  finding  a 
more  practicable  solution  to  the  problem  of  the  induced  velocity  field  about  a 
helicopter  rotor. 
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